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We evaluate the mean value of the van Vleck response of a two-level system with level spacing 
distribution and argue that it describes the orbital magnetism of small conducting particles. 

I. INTRODUCTION 

The significance of the few-level physics for small metal particles has long been realized Q, Q. At temperatures 
smaller than the mean level spacing at the Fermi level, T < A, the spin susceptibility, for instance, is determined 
by the (spin-flip) transitions to the first unoccupied state for the particles that contain even number of electrons 
and is typically much smaller than the Pauli susceptibility. For the odd-electron particles, it is given by the Curie 
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susceptibility. This is in contrast to the T ^> A (and the bulk) case where the leading term is given by the Pauli 
susceptibility ||. Moreover, due to the exponential activation factors associated with transitions to unoccupied states, 
one should expect extremely broad distribution of susceptibility values and, therefore, large variations from particle 
to particle || . 

The orbital magnetic response of small metal particles, on the other hand, does not easily reduce to a few-level 
problem. While it does turn out to be essentially a Fermi level property for T ^> A jj| and in the bulk systems ||- 
j3|, a priori it appears as a property of the entire Fermi sea since all the levels are perturbed by the magnetic field. 
The two competing contributions to orbital magnetism are the precession diamagnetism and the polarization (van 
Vleck) paramagnetism ||. While no rigorous description of the orbital magnetism exists for T < A, it is nonetheless 
I ' expected to be a Fermi level property as well due to the large Fermi sea cancellations between the diamagnetic 
, and paramagnetic contributions. However, at the Fermi level the van Vleck part of the orbital response should be 
dominant and also the one that is very sensitive to variations of the energy level spacings. Consequently, we study a 
model wherein it is assumed that among all the "virtual transitions" between the ground and the excited states of 
the Fermi sea that are responsible for the van Vleck response, the one that determines the particle susceptibility is a 
^} single-electron "transition" from the last occupied to the first unoccupied state. 

Imry has conjectured [0 that, since the electrons in a metal particle should be considered as a canonical ensemble 
JTT| , one can describe the particle response in terms of the difference between the magnetic energies of the canonical 
and grand canonical ensemble Fh — Qh, the latter being the case for a bulk system and relatively well understood. In 
the limiting case of T < A, one should expect the difference to be due to a single electron level crossing the position 
of the chemical potential in the equivalent grand canonical ensemble. However, the formalism based on this approach 
works only in the limit of T 3> A, as explained in Refs. ||, ||, in which case F — = A5N 2 /2 where the particle 
number fluctuation in the equivalent grand canonical ensemble is actually much larger than one, SN 2 ~ T/A 3> 1 
PL In this limit, the perturbation theory can be used p2| . An improvement on this approach involves using the 
exact, non-perturbative level correlation function and, in fact, yields a saturation value of the magnetic energy as T 
approaches A from above. 

In what follows, we first give an argument in support of the large Fermi sea cancellation between the diamagnetic 
and paramagnetic contributions to the magnetic response. We then evaluate the two-level van Vleck response for 
the Fermi level and first unoccupied level. This involves averaging of the inverse energy spacing using the Gaussian 
Orthogonal Ensemble (GOE) statistics for energy eigenvalues and calculation of the magnetic dipolc moment. The 
latter can be estimated using a semiclassical argument but also rigorously evaluated by considering the magnetic 
dipole absorption. We compare our result with the expression for the orbital magnetic response obtained in the limit 
T > A. We restrict our analysis to 2D diffusive particles (disks, narrow strips, and rings) and use the units where 
c = h = 1. We also neglect any sample-specific effects (fluctuations), except when specifically mentioned, and all the 
quantities in consideration are presumed disorder-averaged. 
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II. FERMI SEA CANCELLATION OF VAN VLECK PARAMAGNETISM AND PRECESSION 

DIAMAGNETISM 



The van-Vleck energy of a two-level system in the magnetic field H = Hqz is 
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and the total van Vleck energy is given by 
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where rii = 9 In the case of a continuous spectrum, we can rewrite eq. (^]), using drii/dei = —5 (e^), as 
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where v = (v (0)) is the mean level density at the Fermi level, v (e) is the level density at energy e. We have also used 
the fact that the magnetic moment at the Fermi level can be evaluated semiclassically, that is 
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where r is the classical position vector of an electron and vp is the Fermi velocity. In evaluating the coefficient in eq. 
(||), we averaged over the angle between r and Vf . 

The single-level diamagnetic energy diamagnetic energy is given by 
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and the total diamagnetic energy of the Fermi sea is given by |J 
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It is reasonable to conjecture that the disorder-averaged value (i\ r\ \i) is i -independent which yields, upon converting 
to integration for a continuous spectrum, 



Se 



(tot) 
diarn 



2, ,2 



ve v 



\r\\V)Hl 



16 



(7) 



and is the same as eq. (H), with the opposite sign. This confirms the Fermi sea cancellation between the van Vleck 
paramagnetism and precession diamagnetismQ'Q. However, this derivation does not account for the quantum effects 
beyond the existence of the Fermi sea. Consequently, it is understood that the cancellation is not exact and that a 
Fermi level contribution is not accounted for in the present approximation. The nature of this contribution should 
depend on whether the chemical potential or the number of particles is fixed. In the former case, one expects a Landau 
response, as in bulk systems, while in the latter we make an ansatz of a two-level van Vleck response which involves 
the last occupied (Fermi) level and the first unoccupied level. 

The situation is more complex in a strictly discrete level case where we will give only an order-of-magnitude 

^2 

argument. The first principles estimate of Mif in the diffusive regime is based on the idea first proposed by Shapoval 

jl3| and, later, applied by Gor'kov and Eliashbcrg to small metal particles. Namely, we use the semiclassical 
approach to write 



: This argument can be carried over, with minor modifications, to 3D. 

2 For a disk of radius R, (0| r\ |0) can be evaluated as the area average and equals to R 2 /2. 
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where the bar denotes averaging over all classical trajectories. The correlation time scale r for M (t) is given by 
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where £ = v f t is the electron mean-free-path and D = v 2 F r/2 is the diffusion coefficient. This is because the directions 
of vp are uncorrelated after such time. Also, for Sf ~ Si > r -1 the exponential term becomes oscillatory. The scale 
of r is the relevant sample dimension a and we find, 
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where [Lb is the Bohr magneton and e> is the Fermi energy. Consequently, the order of magnitude value of the 
two-level van Vleck response is obtained from eq. (□) as 
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where \l is the Landau susceptibility |6| and A ~ a 2 is the sample area. The total van Vleck energy can be estimated 
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by multiplying Se^y by (tA) 1 , which determines the range of the integral in eq. (^) , and we find 



Since (0| r\ |0) ~ a 2 , this is in qualitative agreement with eq. (|^). Notice also that eq. (|J) yields 
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for a single-level contribution to the precession diamagnetism. 



III. VAN VLECK PARAMAGNETISM IN THE TWO-LEVEL MODEL 



Turning again to eq. ([!]), where now Si is the energy of the last occupied state (Fermi level) and Ef is the energy 
of the first unoccupied state, we find, averaging with the Wigner-Dyson distribution ]T^ |, 
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Here s is the level degeneracy (s = 2, on the account of spin) and v = sA 1 . We have already estimated the matrix 
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element using the semiclassical argument. However, a precise derivation of Mif in the diffusive regime can be done 

by means of evaluation of the low-frequency magneto-dipole absorption in the field H = Hq exp (— icot) z. In the case of 
the continuous energy spectrum (for instance, when the level broadening 7 is larger than A ), the quantum- mechanical 
expression for the absorption should yield, up to small corrections, the classical value [ fL5| . This is in complete analogy 
with the electric-dipole absorption (barring screening effects for the latter), which is discussed in detail in Ref. p6[ . 

enters into the quantum-mechanical expression (see below), it can be extracted by equating with the 
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classical value of the absorption. 



The classical absorption is readily evaluated according to |17 
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for a disk of radius R and for a narrow metal strip, such that L x <C L y , respectively. Here 
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and a is the Boltzmann conductivity. It is assumed that the frequency is such that 5 3> R,L X . Combining eqs. 
we obtain 
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for the absorption, respectively, in the disk and the strip. 

The quantum absorption, for the continuous spectrum, can be evaluated (in complete analogy with the electric- 
dipole absorption [fl6||) by means of the Fermi golden rule and we find 
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where the small quantum corrections of order A 2 /7 2 (or A 2 /uj 2 if lu > 7) |T^] are neglected. For the electric-dipole 
absorption |l6|] , the quantity corresponding to Mif is the electric dipole matrix element Pjj . The latter can be 
evaluated from first principles in the diffusive approximation and the classical and quantum results can be evaluated 
independently and are equal in the considered limit. On physical grounds, we can assume that this is also true for 
the magnetic-dipole absorption. Consequently, we equate the r.h.s. of eqs. (|l9|), ( p0| ) with that of eq. ( pi] ) and, with 
the use of a = e 2 vD/A, we obtain 
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for the disk and the strip respectively. 



Substituting thus found value of 



if 



in eq. (UJ), we find (for s = 2) the following result for the van Vleck energy: 
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for the disk and the strip respectively. Eqs. (|24j) and (|25|) should be compared with the magnetic part of the energy 
obtained for T > A in the so called "mixed approximation" |] wherein the exact, non-perturbative level correlation 
function is used yet the large particle-number fluctuation in the equivalent grand canonical ensemble is assumed also, 
the latter being true only for T> A, 
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The procedure for the evaluation of r^ 1 using the gauge where the vector potential is tangential to the surface is 
described in Ref. fill and we find 
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for the disk and the strip respectively. 

We also mention the Aharonov-Bohm response of a narrow ring threaded by the flux 0, where the derivation is 
particularly simple. Namely, the van Vleck energy in this case is given by 
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where \vif | 2 is evaluated by considering absorption in the alternating flux <fi exp (~iu>t) which generates the electric field 
according to the Lenz's law that, in turn, produces the Boltzmann current density proportional to v. Consequently, 
we find 

\v,ff = i- v (30) 



and, averaging over the level spacing, 
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The latter has the same parametric dependence as 
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which is the limiting value of energy for T > A || 



IV. DISCUSSION 

The main result of this paper is that the parametric dependence of eqs. (p^), ( p^ ) and (p7|), ( |28| ) is the same 
(with closely matching numerical coefficients) . Consequently, the orbital magnetic response of a level-quantized metal 
particle can be satisfactorily explained in terms of the two-level van Vleck response. This conclusion also holds for 
other level distributions characteristic of disordered (chaotic) systems, such as the Gaussian Unitary Ensemble (only 
the numerical coefficient will be different). 

Whereas the two-level picture might be also valid for the Poisson distribution, which is the case for classically 
integrable systems, the mean response needs to be examined more carefully because level "bunching," characteristic 
of this distribution, implies that for any magnitude of the field there will be a significant number of particles for which 
the perturbation theory no longer applies. 

Another question which we hope to address in a future work is the sample-specific (fluctuation) effects that are 
expected to be quite large in the limit considered here. We point out that, in addition to the fluctuation of the level 
spacing and the magnetic moment, the details of the cancellation between the van Vleck paramagnetism and Landau 
diamagnetism should differ from particle to particle leading to variations in the number of levels contributing to the 
total response. 
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